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Abstract 

We present the asymptotic distribution for two-sided tests based on the profile likelihood 
ratio with lower and upper boundaries on the parameter of interest. This situation is relevant 
for branching ratios and the elements of unitary matrices such as the CKM matrix. 



1 Introduction 

In Ref. [JLJ, the asymptotic distribution of several test statistics based on the profile likelihood ratio 
are presented. In particular, the test statistics for one- and two-sided tests of a single parameter of 
interest /i in the presence of nuisance parameters 6 are given cases where the parameter of interest 
is unbounded and bounded with /U > 0. Here we present the asymptotic distribution for two-sided 
tests based on the profile likelihood ratio with lower and upper boundaries on the parameter of 
interest. This situation is relevant for branching ratios and the elements of unitary matrices such 
as the CKM matrix. 

We consider the case where fi £ Using the notation of Ref. [Ij, the maximum likelihood 

estimate of fi is denoted jl and cr^ = var[/i]. Various approaches to estimating a are presented in 
Ref. [IJ. In order to use the results of Wilks and Wald, the strategy is to consider the situation 
when /i is unbounded and impose the boundary in the test statistic itself. Specifically, the test 
statistic is 

-21nM^ ,-<A<^. (1) 
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where 0{fj,) is the conditional maximum likelihood estimate of 9 given 

2 Result 

The relationship between and fi can be inverted to obtain fi{t^) for the two branches fi < fJ^ and 
fi > fi. Wald's theorem states that 

lexp(-(ii^). ,2) 



Through a straightforward change of variables, one can obtain the distribution for the test statistic 
if^. The pdf f{i^\fj.') is found to be 

fiW) = fL{W) + fRiW) (3) 
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where the dimensionless variables 5- = {fj, — fX-)/a, d'_ = {fi' — fj,-)/a , 5+ = {fi — fj,+)/o 
S'^ = {ji' — lJi+)/o' are used to simphfy the expressions. 

The special case ijl = ijI is therefore 
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6c < h < 

i„ > 6% , 



where 6c = min[// — fi^^jij^ — \j\lo and 6p = max[// — — 
The corresponding cumulative distribution is 
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where $(.x) is the cumulative probability distribution of the standard normal distribution. 
The special case = is therefore 
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3 The critical cutoff 



Confidence intervals are defined as the set of where the test statistic is less than or equal 
to a critical cutoff ka{lj)- The cutoff is chosen to insure the desired coverage probability. For a 
100(1 — a) % confidence level interval, the cutoff is defined by F{ka{^i)\^) = 1 — a. 

When there are no boundaries, the distribution of the test statistic follows a distribution 
and the cutoff is constant. Specifically, it is given by i<^'^il - a/2))^ which gives the familiar 
values of 3.84 for 95%, 2.71 for 90%, and 1 for 68% confidence intervals. 

When will the boundary matter? The critical cutoff is modified for // < //_ +(t<I>^^(1 — a/2) and 
/i > /X4. — cr<I>~^(l — a/2). Thus edges of a confidence interval using the standard cutoff are fine if 
they fall in the intermediate range; however, they will over-cover if they are near the boundaries. 
As a increases the range of /i with a modified cutoff grows. Once a > cJcrit = (/"+ — /^-)/(2\/ia)) 
then there is no region of fi where the cutoff is not affected. 

When testing at the boundary, the critical value kai^i-) = ka{f^+) is always affected. For large 
values oi 6 = — ^-)/a (ie. when fi is well measured with respect to the range of ^) only one 
boundary is important; however, for small values of 5 (ie. when jl is poorly measured with respect 
to the range of jj) both boundaries are important. The cutoff at the boundary is given by 

(2^^^{l-a))5-6'^ (5<<5-i(l-a) 
= ka{n+) = I (11) 
[($-i(l-a))2 <5>$-i(l-a) 

Note that for a 95% confidence interval, if o" > (/i+ — ^_)/1.64, then the far away boundary reduces 
the critical cutoff below the 2.71 one might expect from the presence of the boundary being tested 
and it is significantly smaller than the 3.84 cutoff one has from assuming a distribution neglecting 
any boundary effects. 

Figures [ijjs] show the critical cutoff ka{n) for 68%, 90%, and 95% confidence intervals for a <^ 
(Tcrit, CF = CTcrit, and CJ ~ (/i+ - ;U_)/$~^(1 - a). 



4 Conclusions 

The presence of both lower- and upper-boundaries on a parameter of interest is a common situation 
in particle physics. For example, branching ratios and elements of unitary matrices are bounded 
between and 1. The formulae presented here are essentially the asymptotic versions of the 
Feldman-Cousins approach [2] extended to incorporate nuisance parameters via the use of the 
profile likelihood ratio test statistic as in Ref. [1]. 
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Figure 1: The critical cutoff for 68% confidence level interval for several cases. The upper horizontal 
line corresponds to the naive cutoff neglecting any boundaries of (^•"^(l — a/2))^. The lower 
horizontal line corresponds to the cutoff of ($^^(1 — a))^, which is appropriate when testing on 
the boundary when a is sufficiently small. The dashed curve corresponds to a = (;U+ — /i_)/10 
where the cutoff is only affected for /x < + (T$~^(1 — a/2) and iJ, > iJ,+ — a^~^{l — a/2) (vertical 
lines). The solid curve corresponds to the cutoff for a = (Jcrit; for any value of cr > (Jcrit, the cutoff is 
affected for all values of fi. Finally, the dotted curve shows the case of cr = 1.2{^-f- — fi-)/^^^{l — a), 
where a > — — a) means that the critical cutoff on the boundaries is affected by 

the faraway boundary. 
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Figure 2: The same as Figure [T] for 
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Figure 3: The same as Figure [T] for 95%. 
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